Structure and thermochemistry of K2Rb, KRb2 and K2Rb2 
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The formation and interaction of ultracold polar molecules is a topic of active research. Un- 
derstanding possible reaction paths and molecular combinations requires accurate studies of the 
fragment and product energetics. We have calculated accurate gradient optimized ground state 
structures and zero point corrected atomization energies for the trimers and tetramers formed by 
the reaction of KRb with KRb and corresponding isolated atoms. The K2Rb and KRb2 trimers are 
found to have global minima at the C2v configuration with atomization energies of 6065 and 5931 
cm~^ while the tetramer is found to have two stable planar structures, of D2h and Cs symmetry, 
which have atomization energies of 11131 cm~^ and 11133 cm~^, respectively. We have calculated 
the minimum energy reaction path for the reaction KRb+KRb to K2+Rb2 and found it to be 
barrierless. 



^ ■ 

Oh' 
I 

S 

^ : 

q . 

^ : 

Oh. 



> 

in 

cn 
o 
o 



The formation and interaction of ultracold polar 
molecules is a topic of great current interest in physics. 
New techniques for the formation of rovibrational ground 
state polar molecules via STIRAP[l| (stimulated rapid 
adiabatic passage) or FOPAQ (Feshbach-optimized 
photo-association) allow experiments to be performed 
with w = heteronuclear diatomic molecules, including 
KRb[l-@| and LiCsQ Proposals for quantum computa- 
tion with polar molecules [g, 'o'l have generated a grow- 
ing need for understanding of the dynamics of diatom- 
diatom collisions. Such studies of diatomic dynamics re- 
quire knowledge of the open and closed channels relevant 
in those reactions. The purpose of the present paper is 
to present accurate ah initio calculations of the structure 
and thermochemistry of several chemical species relevant 
to the study of KRb— KRb dimer interactions. 

Theoretical work on electronic structure of few-body 
alkali systems has been limited to lighter homonuclear 
trimers, in particular doublet and quartet [llj Lis, 
doublet K3 [in and quartet Naa [l3|- The recent work 
of Zuchowski and Hutsonjl^ has characterized the atom- 
ization energy of the alkali homo- and heteronuclear tri- 
atomic species formed from Li, Na, K, Rb, and Cs. These 
homonuclear trimers have A' ground electronic states in 
Cs symmetry that correlate to B2 symmetry in C2V ■ Pre- 
vious mixed alkali tetramer studies have been limited 
to structure studies of Li„Xm (X=Na and K) [l^, [T6| 
and that of RbCs-|-RbCs To date no such calcula- 
tions have been reported for the heteronuclear K„Rbm 
tetramer molecules. 

Electronic structure calculations were performed on 
K2, Rba, KRb, KsRb, KRbz, and KaRba at the 
CCSD(T) [18] level of theory. As core- valence effects 
can be important in alkali metals, we correlate the inner 
valence electrons in potassium, keeping only ls^2s^2p^ in 
the core. Rubidium is heavy enough that relativistic ef- 
fects are significant, so we replace its inner shell electrons 
by the Stuttgart small-core relativistic (ECP28MDF) 
ECP [l^ . Basis sets are taken from the Karlsruhe def2- 
TZVPP and def2-QZVPP [21^ orbital and fitting 
sets. 

Optimized geometries for K2, Rba, KRb, KaRb, KRba, 
and KaRbs were found at the CCSD(T)/def2-TZVPP 



level of theory. Calculation of the harmonic vibrational 
frequencies was done to verify that the calculated struc- 
tures were minima on the potential energy surface, and 
the calculated frequencies were used to obtain vibrational 
zero point energy (ZPE) corrections. These structures 
were further optimized at the CCSD(T)/def2-QZVPP 
level of theory, leading to a 0.07 A correction in the 
bond lengths and 60 cm^^ in final atomization energies. 
The CCSD(T)/def2-QZVPP geometries are tabulated in 
Table HI 

Evaluation of the contribution of scalar relativistic cor- 
rections to K2 indicate a small 0.005 A and < 8 cm~^ 
contribution in all electron correlation calculations [24i]. 
while for Rba it has been shown 25] that the small core 
Stuttgart pseudopotential gives an accurate representa- 
tion of relativistic effects on the bond length and disso- 
ciation energy. 

Single point energy calculations were then done us- 
ing the CCSD(T)-F12b [H, ^ (explicitly correlated 
CCSD(T)) level of theory. The use of explicitly corre- 
lated methods accelerate the slow convergence of the one- 
particle basis set by including terms containing the inter- 
electron coordinates into the wavefunction [28*], thus 
yielding very accurate results using triple and quadruple 
zeta basis sets. In addition, we estimate the complete 
basis set (CBS) limit usiiig the two-point extrapolation 
formula of Helgaker et al [23] 



rv^E^ - (n - ifE^.i 
rfi — {n — 1)^ 



(1) 



In Table ITT] the CCSD(T) and CCSD(T)-F12b dissocia- 
tion energies for the def2-TZVPP and def2-QZVPP basis 
sets are tabulated as well as the zero point energy (ZPE) 
corrected atomization energies. After extrapolation, the 
diatomic CCSD(T)-F12b ZPE corrected dissociation en- 
ergies agree very well with the experimental diatomic dis- 
sociation energies, as shown in Table [III The ab initio 
calculations were done using the Gaussian 09 [s^] and 
MOLPRO [31U33I packages. 

We have found that both KaRb and KRba have two 
energetically close local minima on the ground state sur- 
face, one of Cat; symmetry and another less symmetric Cs 
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TABLE I. Calculated CCSD(T)/QZVPP molecular geometries (in Angstroms and degrees). 



Ka 3.956 
Rba 4.233 
KRb 4.160 





?'K-Rb 


'"K-Rb 


e 




KaRb C2„ 


4.279 


4.279 


70.68 




KaRb Cs 


4.361 


5.234 


48.81 




KRba C2„ 


4.271 


4.271 


82.13 




KRba Cs 


4.193 


5.179 


57.07 






»'Rb-Rb 




^K-Rb 


fe-Rb-Rb fe-K-Rb 


KaRba D^h 


8.224 


4.0307 


4.579 




KaRba Cs 


4.761 


4.408 


4.189 


53.34 55.476 



TABLE IL Dissociation and zero point energies calculated using CCSD(T) and CCSD(T)-F12b correlation methods with 
successive basis sets and CBS extrapolated values (in cm~^). 





ZPE TZVPP 


De TZVPP 


De QZVPP 


Do 


CBS 




CCSD(T) 


CCSD(T) 


CCSD(T)-F12b 


CCSD(T) 


CCSD(T)-F12b 


CCSD(T) 


CCSD(T)-F12b 




46.0 


4098.8 


4276.9 


4460.0 


4369.7 


4677.6 


4391.5 


Rb:^ 


26.8 


3494.3 


3723.3 


3842.7 


3885.4 


4070.2 


3976.8 


KRb'^ 


35.4 


3829.4 


4015.6 


4135.6 


4128.7 


4323.6 


4175.7 


KaRb C2„ 


69.8 


5588.2 


5805.5 


6067.7 


5995.7 


6574.2 


6009.4 


KaRb Cs 


72.4 


5606.3 


5843.7 


6179.1 


6015.9 


6524.7 


6069.1 


KRba Ca„ 


62.8 


5394.5 


5635.1 


5911.0 


5842.2 


6043.5 


5788.3 


KRba Cs 


59.0 


5215.9 


5475.4 


5728.5 


5690.4 


6225.1 


5930.5 


KaRba D^h 


129.5 


10210.8 


10669.4 


11275.3 


11011.1 


11922.7 


11131.0 


KaRba Cs 


126.2 


10198.3 


10629.9 


11211.4 


10993.7 


11824.6 


11133.0 



Experimental value 4405.389 cm ~^22i] . 
^ Experimental value 3965.8 cm~^ [2a[. 
Experimental value 4180.417 cm~^[3|. 



structure (geometries given in Table H]). While dependent 
on the level of theory used to evaluate the atomization 
energy, we conclude that the symmetric C^v geometry is 
the global minima for each trimer. The atomization en- 
ergies calculated are found to be in good agreement with 
those recently published by Zuchowski and Hutson 14 1. 

The K2Rb2 tetramer is found to have two nearly de- 
generate minima on the potential energy surface. One is 
a rhombic structure of Z?2/i symmetry, and another pla- 
nar (Cs) structure that corresponds to an interchange of 
K and Rb atoms. These structures are bound by ~ 3000 
cm~^ with respect to K2-|-Rb2 or KRb-|-KRb. The elec- 
tronic structure of these two isomers is very similar, and 
their stability is likely due to three-center bonds of the 
sort proposed for Li„Na4_„ clusters [TsUlGj . The rhombic 
K2Rb2 structure has a short (~ 4A) distance and a long 
(^ SA) Rb-Rb distance. The equivalent structure where 
the K-K distance is short and the Rb-Rb distance is long 
is found to be a transition state, not a stable minimum. 

To determine if there is any barrier to the 
KRb-|-KRb^-K2Rb2 — >Rb2-l-K2 reaction, we calculate a 
minimum energy path for the KRb-|-KRb— ?'K2Rb2 and 
Rb2-f K2— J>K2Rb2 reactions. We start by locating the 
minimum energy geometric configuration at long range. 



This is done by calculating ah initio the dipole and 
quadrupole electrostatic moments of K2, Rb2 and KRb 
and then minimizing the long range electrostatic interac- 
tion energy [34] with respect to the angular configuration 
of the molecules. This minimization resulted in a T type 
geometry for both K2-I-Rb2 and KRb-|-KRb. We have 
recently shown that long-range expansions of this type 
accurately reproduce diatom-diatom interaction energies 
[35| . From these initial geometries, the reaction path 
was followed by freezing the diatom-diatom distance and 
optimizing the diatomic bond lengths and angular orien- 
tations at the frozen core CCSD(T)/def2-TZVPP level 
of theory. Single point energies were evaluated along this 
path using the CCSD(T)-F12b level of theory including 
the core-valence correlation energy and extrapolated to 
the CBS limit as discussed above. This procedure, in 
which a high level energy profile is evaluated along a re- 
action path calculated at a lower level of theory, is known 
to be a good approximation to the energy profile along 
the reaction path calculated at the high level of theory 

Hi. 

We find that the KRb+KRb dissociation limit con- 
nects to the minima while the K2-I-Rb2 dissociation 
limit connects to the Cg minima, with no barrier found 
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FIG. 1. Minimum energy path connecting the KRb-|-KRb 
and K2+Rb2 dissociation limits. Included are schematic ge- 
ometric at points of interest, where open and closed circles 
represent rubidium and potassium atoms respectively. 



to either reaction. A similar conclusion was obtained 
for the RbCs-f Rbcs— )>Rb2-l-Cs2 reaction by Tscherbul et 
aZ(T^. To finish characterizing the reaction path go- 
ing from dissociation limit to the other, we locate the 
transition state and calculate the intrinsic reaction coor- 
dinate (IRC) L37| reaction path connecting the Cs and 
D2h minima structures at the same level of theory as de- 
scribe above. Optimizing the transition state geometry at 
the inner valence CCSD(T)/def2-TZVPP discussed pre- 
viously and evaluating an accurate atomization energy 
using our CCSD(T)-fl2b prescription we find that the 
transition state is 1167.3 cm~^ above the dissocia- 
tion energy. The calculated reaction path is plotted in 
Figure [1] using the approximate reaction coordinate 

AR — (i?R,b-Rb + ^Cs-Cs)/2-(i?Rb-Cs + -RRb-Cs)/2 (2) 

where Ra-b is the distance between atoms A and B. 

The formation and trapping of rovibrational ground 
state KRb diatoms with a high phase space densityjl] 
offers the opportunity to study chemical reactions in the 
ultra-cold regime [6]. As seen in Figure the three-body 
reaction KRb+Rb— ^Rb2+K is energetically forbidden at 
ultra-cold temperatures, leaving the endothermic four- 
body reaction KRb-f KRb^-Rb2-f K2 as the only path- 
way to forming Rb2 within the trap. Measurements of 



the population of Rb2 within the trap will then allow di- 
rect probing of the exchange reaction rate of KRb-|-KRb. 
Inherent in this exchange reaction is the possibility of 
studying the role of fermionic/bosonic spin statistics in 
ultra-cold chemical reactions ,38i - l43 | . In this tempera- 
ture regime, s-wave scattering of fermionic ^"^KRb is sup- 
pressed which should greatly diminish the reaction rate 
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FIG. 2. Schematic energy level diagram for fragment and 
structure energies involving KRb with KRb and separated 
atoms. Inset figure shows the small difference between the 
KRb-fKRb and K2+Rb2 asymptotes. 



of ^°KRb-t-'*'^KRb, thus leaving the trap stable to four- 
body losses. If instead the trap was formed with bosonic 
■^^KRb or ^^KRb molecules, no such coUisional suppres- 
sion is expected, where we then expect comparably large 
reaction rates to occur. It is also possible to explore 
recent theoretical predictions f44] which show that if a 
bosonic dimer is composed of two fermions of very differ- 
ent masses the resulting exchange reaction should still be 
suppressed despite the overall bosonic nature. This could 
be accomplished by using fermionic ^''K and a long lived 
®*Rb or ^^Rb. The comparison between reaction rates 
in the above described interactions can then be used to 
directly study the effects of fermion/boson spin statistics 
to that of chemical reactions. 
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